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Abstract 

Analysis of semi-inclusive DIS hadroproduction suggests broadening of transverse momentum 



Q-i! distributions at x below a few 10 ^ , which can be modeled in the Collins-Soper-Sterman formalism 

Q I by a modification of impact-parameter-dependent parton densities. We discuss the consequences 



of such a modification for the production of electroweak bosons at hadron-hadron colliders. If 



^ . substantial small-x broadening is observed in forward boson production in the Tevatron Run-2, 

\ it will strongly affect predicted qx distributions for W^, Z^, and Higgs boson production at the 

Large Hadron Collider. 
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I. INTRODUCTION 
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Moderately heavy particles with masses about 100 GeV can be copiously produced in 
hadron collisions at TeV energies. Production and observation of such particles is accompa- 
nied by intense hadronic activity, which affects overall rates and differential cross sections. 
Radiation of several hadrons of relatively low energy may produce a strong recoil effect on the 
heavy particle. Calculations in perturbative quantum chromodynamics (PQCD) may have 
to be reorganized to resum large logarithms associated with multiple hadronic radiation. 

Electroweak symmetry breaking in the standard model introduces vector bosons and 
Z'^ with masses of about 80 and 91 GeV, and a scalar Higgs boson with a likely mass below 
250 GeV, as suggested by global fits of electroweak data. Supersymmetric extensions of the 
standard model require existence of several types of scalar bosons, with the lightest naturaU 
having a mass below about 140 GeV in the minimal supersymmetric standard model 
Identification of heavy bosons relies on understanding of transverse momentum distributions 
for heavy bosons and their decay products, as well as for background processes. Short-lived 
bosons are commonly detected by observing their decay into a pair of color-neutral particles 
with large transverse momenta (of order of a half of the boson's mass) in the laboratory 
reference frame. Furthermore, the mass My of the heavy boson (e.g., W boson) can be 
determined from the transverse momentum distribution of the decay products. 

When the boson's transverse momentum Qt is much smaller than the boson's virtuality 
Q, the calculation of the transverse momentum distribution must evaluate an all-order sum 
of large logarithms In"" (qT/Q)- The resummation of small Qt logarithms can be realized using 

onnnnQ 

one of several available methods; see [3, y, U, I^UO, lZ| for some recent references. As stated 
by Collins, Soper, and Sterman (CSS) in Ref. (sj, at moderate energies all large logarithms 
are generated by a product of a soft (Sudakov) exponential and special parton distribution 
functions in the impact parameter (6) space. The 6-space resummed form factor can be 
evaluated in PQCD at large Q {Q"^ ^ ^qcd) small b { <^ Aq^^). Its calculation 
is possible due to the collinear factorization of hadronic cross sections, valid when Q is 
not small compared to the total energy of the hadronic collision. The CSS approach 
describes well the available qt data from fixed-target Drell-Yan experiments and Z boson 
production at the Tevatron 
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In this paper, ^ we discuss possible deviations from this established picture that may 
occur in the new range of energies accessible to the Large Hadron Collider (LHC). As we 
move from the 1.96 TeV pp collider Tevatron to the 14 TeV pp collider LHC, the typical 
fraction of the collision energy going into production of moderately heavy bosons decreases 
below one percent, leading to several new effects. Transverse momentum distributions of 
W and Z bosons can be altered at small momentum fractions x by increased contributions 
from qg and gg hard scattering, which tends to produce electroweak bosons with larger 
transverse momenta as compared to the dominant process of qq scattering. One also expects 
enhancement of the logarithms ln(l/x) in the matrix elements of order and beyond, the 
effect not completely accounted for in the conventional qx resummation. Such logarithms 
may be enhanced by a larger QCD coupling strength as at qr less than a few GeV, even 
when they are negligible in inclusive cross sections depending on one hard scale Q ^ I GeV. 
Finally, complex nonperturbative dynamics contributing at 6 > 1 GeV~^ may also depend 
on X. At present, the magnitude of the x-dependent corrections to W, Z, and Higgs boson 
production at the LHC energy is largely unknown, in part because no experimental data on 
qr distributions is available yet in Drell-Yan-like processes at x of a few 10"'^ or less. 

Resummation of qx logarithms at qr/Q and x is realized in the Dokshitzer- 
G„bov.L,atov.AUa.m.Pa..i (DGLAP) ,c..e of had.on. s.aUenng UUUU, 
which assumes that the dominant parton emissions are strongly ordered in transverse mo- 
menta. At asymptotically high collision energies (x ~ Q/y/S with qr/Q fixed), the 
collinear DGLAP factorization is superseded bv the off-shell factorization in the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) framework jl^, 12, 3|, which sums up leading logarithms of 
1/x, while including g-r logarithms at a finite order of only. When neither DGLAP nor 
BFKL dynamics dominates (for instance, when both x and qx are small), methods for simul- 
taneous summation of qr and 1/x logarithms (e.g., the Ciafaloni-Catani-Fiorani-Marchesini 
e,.a«o„ mm be needed. Together with .be increased .a,„i.,de of « and 
gg scattering, reduced convergence of the perturbation series in the DGLAP framework may 
signal transition to the small-x dynamical regime at the existing colliders. 

Some of the transition phenomena may be already observed in part of phase space at 

^ Our results were partly presented in Ref. 0| as a contribution to the 3d Les Houches workshop "Physics 
at TeV colliders" [JJ]. 
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the ep collider HERA. The analog of the Drell-Yan process in lepton-nucleon deep inelastic 
scattering is semi-inclusive production of hadronic final states, e + p e + hadron(s) + X. 
The CSS resummation formalism can be applied to semi-inclusive deep inelastic scatter- 
ing fSIDIS) to describe the dependence on the Lorentz-invariant transverse momentum qt 

n ~i nn 

|23l . 12J|. Refs. |2^, l2a| have compared predictions of qt resummation to the data for the 
transverse energy flow in the current fragmentation region l^^- The experimentally 
observed qt distribution at x below 10"^ becomes wider as x decreases. This "broadening" 
of qt distributions was modeled by including an extra x-dependent term in the resummed 
form factor. To describe the data, the extra term must grow quickly as x — 0. It provides 
a phenomenological parametrization for all substantial factors responsible for x dependence 
beyond that included in the 0{as) resummed qt cross section. For example, large correc- 
tions from hard scattering of order O(a^ ) substantially modify the SIDIS cross section at 
large transverse momenta [3, 5 corrections are also likely important at 

small qr (not discussed in the published studies). Apart from the hard-scattering contri- 
butions, the dynamics of collinear radiation may change at moderately small x {x < 10^^) 
and small energy scales (1/& < 1 — 2 GeV) as a result of preasymptotic suppression of gluon 
splittings by BFKL logarithms Both enhancement of perturbative scattering at larger 
qx and preasymptotic BFKL suppression of collinear radiation at smaller g-r would result in 
broader distributions. 

While a systematic framework describing all x-dependent effects is yet to be developed, 
some predictions can be made by exploiting universality of soft and collinear radiative cor- 
rections, which dominate the SIDIS energy flow data at < 2 — 3 GeV. For this reason, we 
employ the phenomenological parametrization for the small-g-r cross section found in SIDIS 
to predict the x dependence of g-r distributions at x < 10~^ at hadron-hadron colliders. The 
proposed new form of the Drell-Yan resummed form factor at x < 10"^ is given by 

W{b, Q) = WsLNYib, Q)e-P^-A)b^-Pi-B)y\ (1) 

where WBLNvib, Q) is the resummed form factor found in the global fit 0| to the Drell-Yan 
data at larger x, and the exponential Q~p'y^A)h^-p(xB)b'^ parametrizes the small-x broadening. 

New phenomena at small momentum fractions may have consequences for precision mea- 
surements and searches for new physics. Modifications in the transverse momentum distri- 
butions for W and Z bosons will affect the measurements of the W boson mass and width, as 
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well as the W and Z boson background in the search for new gauge bosons, supersymmetry, 
particle compositeness, extra spatial dimensions, etc. At the LHC, the qr broadening may 
affect detection of the Higgs bosons by altering qx distributions of Higgs bosons and the 
relevant QCD background. 

We examine possible small-x effects in W, Z, and Higgs boson production at the Tevatron 
and LHC. Sizable qr broadening may occur at both colliders. The validity of our model can 
be tested in the immediate future by the analysis of forward rapidity qr distributions in the 
Tevatron Run-2. An observation of signs of qr broadening at the Tevatron will suggest a 
large magnitude of this effect at the LHC. 

In Section HIl we propose a relationship between the SIDIS and Drell-Yan qr distribu- 
tions based on factorization properties of the CSS resummed cross section. To describe the 
enhanced x-dependent contributions, we include in the resummed form factor an additional 
phenomenological term determined from the SIDIS data. In Section IIIIl we predict within 
this model the qr distributions in production of W and Z bosons in the forward region at 
the Tevatron and the central region at the LHC. We also explore implications for light Higgs 
boson production via gluon-gluon fusion at the LHC. The appendix reviews the rapidity 
coverage at the existing hadron colliders. 



II. RELATIONSHIP BETWEEN THE RESUMMED CROSS SECTIONS IN 
DRELL-YAN AND SEMI-INCLUSIVE DIS PROCESSES 

A. Resummed form factor in electroweak boson production 

We start by discussing production of the electroweak bosons at moderate collision energies 
v^, when the ratio Q/VS is large, but not very close to unity. At such energies, the 
logarithms ln(x) or ln(l — x) of the partonic momentum fractions x are of a moderate 
magnitude in most of the events. This situation is typical for the production of Drell-Yan 
lepton pairs in fixed-target experiments and W or Z bosons at the Tevatron. 

At small transverse momenta qr, the differential cross section takes the form 

/* OO ^^^^ 

/ — Jo{qTb)W{b,Q,XA,XB) + Y{qT,Q,XA,XB), (2) 
Jo 27r 



dQ^dydqj 



iT 

where y = (1/2) In [{E + pz) /{E — p^)] is the rapidity of the vector boson, xa,b = Qe^^ / 
are the Born-level partonic momentum fractions, and Joiqrb) is the Bessel function. The 



integral is the Fourier-Bessel transform of a form factor W{b,Q,XA,XB), given in impact 
parameter (b) space. It contains an all-order sum of the logarithms a"ln'"(g|,/(5^), which 
dominates in the small-gr region. Y{qT,Q,XA,XB) is the regular part, defined as the dif- 
ference of the fixed-order cross section and expansion of the Fourier-Bessel integral to the 
same order of a^. It is numerically small at — * 0. 

W{b,Q,XA,XB) can be expressed as a product of Born-level prefactors cr^°'*, a Sudakov 
exponent e"'^^''''^^ and 6-space parton distributions P„ {x,b): 

Wib, Q, XA, Xb) = ^Y1 ^''^'''^^ 'Pa\xA, b) Vf\xB. b) . (3) 

a, 6 

This form will be referred to as "the CSS representation". It has been stated as a factorization 
theorem in Ref. and recently proved in Ref. The summation in Eq. (jH]) is over the 

relevant parton flavors a and b. In the approximation of a vanishing width of the boson, the 
non-zero Born-level prefactors are given by 

af = ^V2M^Gp \V.,f 6{Q' - M^) (4) 
in W boson production for i = u,d, j = u,d,...; 



(0) 



^V2MlGp ((1 - 4 |e,| sin^ 6^)' + l) SiQ'' - M|) (5) 



in Z boson production for i = u,d, and 



a. 



■V2Gf5{Q^ - Ml) (6) 



(0) 

5767r 

in Higgs boson production. Here M\y, Mz, and are the masses of the VF, Z, and Higgs 
bosons, Gp is the Fermi constant, 6yj = sin^^ (1 — M^/M|) is the weak mixing angle, 
Vij is the Cabibbo-Kobayshi-Maskawa matrix, and Cj = 2/3 or —1/3 are the fractional 
quark charges for up- or down-type quarks. The prefactor (P) for Higgs boson production 
is derived bv using the effective ggH vertex in the approximation of inflnitely heavy top 
quarks This approximation works well for the Higgs bosons that are lighter 



than the top quark. 



^ The factorization for the resummed form factors in the Drell-Yan and SIDIS processes has been also 
demonstrated recently within the framework based on gauge-invariant 6-dependent parton distributions 
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The Sudakov term S{b, Q) is an integral of the functions A {as{fi)) and B (Q;s(/i)), which 
can be calculated in PQCD when b is small (6^ <^ Ag^^): 



S{b,Q) 



2 



Q 



2 



A In — + i3 {asifi)) 



(7) 



bo = 2e~'^^ is a constant parameter involving the Euler constant, ■je = 0.577... . Approxi- 
mation of 5(6, Q) at a finite order of as yields "the perturbative Sudakov factor" Sp{b,Q). 
The perturbative coefficients A^^^ and i?^*-* for the functions A (ajj2)) and B {as{j2)) up to 
the next-to-next-to-leading order are published in Refs. 0, 140I for the Drell-Yan process, 
and Refs. jlll, Q Q, Q for Higgs boson production.^ 

The form factor (pj depends on the light-cone momentum fractions xa and Xb through 
the functions {x,b). In the CSS representation, we distinguish between the functions 

(x, 6) and {x, b) appearing in reactions with timelike (T) and spacelike {S) elec- 
troweak bosons. The relationship between (a;, b) (relevant to the Drell-Yan process) and 

{x, b) (relevant to SIDIS) is discussed later in this section. Expansion of {x-, b) in 
leading powers of 6 at 6 —> gives 

vTi;-M = E/' j^T' (f.M, (^)) /<: (f . ^) + ... - [ci-f ® /.] fet) + ... .(8) 

Here C^^jf"' {x,b,as{bQ/b)) are the Wilson coefficient functions for incoming ("in") partons, 
and fc{x, bo/b) are the conventional parton distributions (integrated over the parton's trans- 
verse momentum k-r). The ellipses denote the terms that are suppressed by positive powers 
of b. The renormalization and factorization scales on the right-hand side of Eq. ^ are equal 
to bo/b. Whenever the series in powers of as{bo/b) and b in Eq. (jH)) converge, ix,b) can 
be approximated by neglecting the positive powers of 6, evaluating C^^jf^ at a finite order of 
cts, and using a parametrization for fc{x,^) from a phenomenological global fit. Estimates 



)m a pnenc 



for C^i'^^^ of order C(as) can be found in Refs. 

Contributions from nonperturbative impact parameters (6 > 1 GeV~^) are suppressed 
at Q 00 by the shape of 14^(6, Q, x^, x^) and oscillations of Jo^qrb). Nonetheless, mild 
sensitivity to the large-6 behavior remains in W and Z boson production when qt is of order 
a few GeV. The factorization of W{b, Q, xa, xb) predicts universality of the nonperturbative 
terms within classes of similar processes, such as production of Drell-Yan pairs, W, and Z 



Equations in this paper correspond to the "canonical" choice |3| of the momentum scales 



bosons. While the nonperturbative terms cannot be calculated yet from the first principles, 
they can be introduced in the calculation by using one of the available models. The present 
Drell-Yan and Z boson data at x above a few 10~^ is consistent with universality, as it was 
demonstrated recently within the "6*" model 0,131 for the nonperturbative terms. The 
"6*" model parametrizes the resummed form factor as 



W{h, Q, XA, xb) = W^'''-\K, Q, XA, xb) e 



-SNp{b,Q-X) 



(9) 



where 



W^''-\K,Q,xa,xb)^^Y1 



(0) 
a,b 



a 



xa, 



(10) 



is the perturbative approximation for 1^(6, Q, xa, x^), evaluated as a function of 6* = 6(1 + 
^'^ /^max)~^^'^ for bmax ~ 1 GeV~^ The nonperturbative Sudakov function, 



SNp{b,Q]b^) = - In 



W{b,Q,XA,XB) 



(11) 



WP^^%b,,Q,XA,XB) 

is determined from the experimental data. If the flavor dependence of large-6 contributions 
is neglected (a reasonable first approximation in reactions dominated by the scattering of u 
and d quarks), Sj\fp{b, Q; 6*) is straightforwardly related to the true Sudakov factor S{b, Q) 
and 6-dependent parton distributions ixA,b) in Eq. (jH]): 

SNp{b, Q; &.) ^ S{b, Q) - Sp{b,, Q)-\n\ ^l! ] - In ^ ^^"^^'^^ 



SPf\xB, b^] 



(12) 



This representation realizes the generic form for SNp{b, Q; 6^,) found in Ref. 

SNp{b,Q]b^) = 5(5(6, In Q; 6,) +ga{b,XA;b^) + gb{b,XB]b^), 



where gs, ga, and gi, are some functions of the shown arguments. Several analyses have 
)een carried out in the past in order to constrain SMp(b,Q:b^) from the Drell-Yan data 
2mm. the p.se. wo.k, we the BLNY p...e.i..ion Q, given 
by 

Q 



qBLNY 
'->NP 



5, Q; 



0.21 + 0.68 In 



0.126 \n{100xAXB) 



(13) 



3.2 GeV^ 

for bmax = 0.5 GeV~^. This choice of parameters yields good agreement with the Drell-Yan 
data and is convenient for comparison with the qt fit in SIDIS, which was also performed 
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for bmax = 0.5 GeV ^. Alternative forms for the parametrization ofW{b,Q,XA,XB) at large 



Q,y,yi. 



b have been also considered Studies of the general structure of the large- 

b contributions with the methods of infrared renormalon analysis have demonstrated, in 
particular, that the leading power-suppressed contribution is quadratic in b (Gaussian), i.e., 
that Snp ~ 6^, as in Eq. 



Since xaXb = Q'^/S, S 



BLNY I 
NP \ 



'j,Q;b^:) is independent of the boson's rapidity. It will be 



used in Section UTTl to predict the W and Z boson resummed cross sections in the case when 



all rapidity dependence is contributed by 



C 



in (T) 
a/c 



fc 



(x,6o/&*)- By comparing Eqs. (|T^ 



and (fT3|l . we find that the logarithmic combination of (x, 6), 



In 



In 



V?{x,b) 



C 



in (T) 



X 



bo 
' b. 



a/c /c 

is a slowly varying function of x (denoted by c^'^\x, b; b^)). Therefore, 

^^ cm(x,b-M 



(14) 



pin (T) 
'-a/c 



fc 



(15) 



when X is larger than a few percent. 

The nonperturbative function SNp{b, Q; 6*) may take a different form in gluon-dominated 
Higgs boson production. Fortunately, the sensitivity to SNp{b, Q; 6*) is suppressed in gluon- 
initiated channels at the LHC by a larger color factor associated with the perturbative 
leading-logarithm contributions. In Section UlI C 2\ we quantitatively estimate the effect of 
the uncertainty in S^pib, Q; 6*) on the Higgs boson cross section. 



B. Resummed form factor at small x 

As X 0, the approximation (fT5| for "P^ (x, b) may be altered because of the re- 
duced convergence of the perturbation series at intermediate and large impact parame- 
ters. The factorized expression Q for W{b,Q,XA,XB) is derived from the requirements of 
renormalization- and gauge-group invariance, and does not rely on the perturbative conver- 
gence for the individual components, including (x, b). Similar factorization holds for the 
resummed form factor in semi-inclusive deep inelastic scattering (SIDIS), e+A e+B+X 



Q, y, y, 



where A and B are the initial- and final-state hadrons, respectively. 
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SIDIS probes the x-dependence of the spacelike distributions Vj {x, h) for quark flavors 
(j = M, d, ...), which differ from the timelike quark distributions Vj (x, h) by a factor r{h): 

Vf\x,b) =r{b)vf\x,b). (16) 



This relationship is derived in Section Hi D| which also shows that r(b) is independent of x. 
The experimentally observed hadronic energy flow in SIDIS l^,!^! is compatible with 



the following form of vf\x, b) (2^, [2^ at 10"^ < a; < 1: 



3 



10-^<x<l: Vf\x,b) 



pin (5) 



fa {x 6*) e"°-°^^''^/''+''*^'(^''''^*) 



(17) 



1 



-C 



(18) 



with 

C'"^'^Vt b 1 - ^ 

SVa ^'^'''*^-r(6,)^> 

The exponential exp (— 0.0136^/x) results in the additional x-dependent broadening of the 
qr distributions. The power of the exponential, 0.0136^/x, is small at x 10"^ and rapidly 
grows at 10"^ < x < lO'^. The SIDIS flt does not constrain vf\x, 6) at x < 10"^. The 1/x 
growth of the exponential power may be modifled at very small x by the turn-on of other 
scattering mechanisms, notably BFKL resummation and saturation. The remaining terms 
(collectively denoted as c^^\x^b]b^)) do not depend strongly on x and can be numerically 
important only at large x and b. The behavior oiV j (x, b) in the transition region x ~ 10~^ 
and precise form of c'^^\x,b]b^) are not well determined by the SIDIS fit because of the 
limited accuracy of the data. 

According to Eq. (fTHll . the broadening exponential must also be present in the timelike 
distributions Vj (x, b) and dominate at moderately small x: 

J J (^^^^j g-o.oi3feV^+c(^H^',6;b.)_ 

The functions c^'^\x, b; 6*) and c^^\x, b; b^) are in principle related by 

c(^)(x,6; b,) = c^^\x,b- 6*)+ In 



10"' <x < 10-2 : Vf\x,b) 



pin (T) 



r{b) 



(19) 



(20) 



r{b*)\ ■ 

We find that the large-x form ^ of the resummed form factor is quantitatively known from 
the fits to the Drell-Yan and Tevatron Run-1 Z boson data, while the leading behavior at x < 
10~^ can be inferred from SIDIS. What is not precisely known is the behavior of V j (x, b) 
in the transition region around x ~ 10~^. To estimate the possible impact on the Tevatron 
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and LHC observables, we interpolate Vj (x, b) between the large- a; parametrization (|T 
and small- a; parametrization (fTHll by using a smooth trial function p{x): 



10-^<x< 1 : Vf\x,b) 



Mn (T) 



fa 



(21) 



We choose p(x) as 

p{x) = cA\ -^ + \--\, (22) 

y y x5 xqJ 

where Cq controls the magnitude of the broadening for a given x, and Xq is a characteristic 
value of X below which p{x) becomes non-negligible. In the limits x ^ xq and a; <^ xq, 
Eq. (OT|l reduces to the large-x form (fTK| and small-x form (fTn|l , respectively. As a result of 
the new form for vf"\x,b), the BLNY representation Q is multiplied by ^-p('-^A)b^-p{xB)b'^ . 



W{b,Q,XA,XB) 



a,b 



pin (T) 
'-a/c 



fc 



bo 

0* 



r,in (T) 



fd 



bo 

xb, 7- 

0* 



X e 



-Sp{b,,Q)-SBL^^{b,Q;K) -b^p{xA)-b^p{xB) 



(23) 



Eq. ^ is restored from Eq. when both x^ and Xb are much larger than xq. The 

Gaussian exponential e"''^^')*^ approximately preserves the integral of the Qt distribution, 
and the observables inclusive in qx remain essentially unaffected. To reproduce the leading 
small-x behavior in Eq. (fT9|) . our quantitative estimates for W and Z boson cross sections 
in Section IIIII are made for Cq = 0.013 and Xq = 0.005. We will also consider variations 
around these values to examine sensitivity of the Tevatron W and Z production to various 
broadening scenarios. The cross sections for Higgs boson production are computed for two 
values of cq, 0.013 and 0.026, in order to evaluate the uncertainty due to the unknown small-x 
form of the gluon distribution V„ (x, 6) (not constrained by the SIDIS data). 



C. Resummed z flow 



7* (9) 

In the remainder of this section, we review the resummation in SIDIS e{k) + A{pa) — > 
e{k') + B{pb) + X and its relationship to the resummation in the Drell-Yan process. The x 
dependence of the distribution V„ (x, b) can be determined from the data on the hadronic 
2;-flow Tjz, defined by 3, | 



= ^ /" z a{e + A ^ e + B + X) dz. 



(24) 
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The z-flow is obtained from the product of the cross section and final-state light-cone variable 
^ = {Pb ■ Va) /{q ■ Va) by integrating over z and summing over all final-state hadrons B. It 
is related to the pseudorapidity (77) distribution of the transverse energy flow {Et) in the 
center-of-mass frame (cm.) of the virtual photon and initial-state hadron A (with 7* moving 
in the +z direction): 

dxdQ^dqT dxdQ'^drj^ 
where = —q^ > 0, and x = Q'^/{2pA ■ q)- The Lorentz-invariant transverse momentum 
qT is a function of rj, qx = Q{x^^ — l^^'^e'^ The resummation applies in the limit 

qT — > 0, which corresponds to the current fragmentation region in the photon-hadron cm. 
frame {rj +00). In this limit, the enhanced soft and collinear logarithms in the theory 
prediction must be summed to all orders. 

Similarly to the Drell-Yan case of Eq. Q , the resummed 2;-flow is given by a combination 
of the Fourier-Bessel transform integral and finite term Yz'. 

poo t.Jt. 

= 2 -j„MWAkQ,.) + y.(,M (26) 

Under HERA conditions, the finite piece Yz{qT, Q, x) is small at < 2 — 3 GeV. The analysis 
of the SIDIS z-flow at such qx probes universal soft and collinear radiation described by the 
resummed Fourier-Bessel integral. The form factor Wz{b,Q,x) is structurally close to the 
resummed form factor (jH]) in the Drell-Yan process: 

Wz{b, Q,x) = ^ Yl e-^(^'«)cr*(6)pf (x, b), (27) 

j=u,u,d,d,... 

where SeA = (k + pA)^, 

af . (1 - ^ + ^) el (28) 

SeAX^ V VDIS VdIsJ 

and i/Dis = Q"^ I {xSeA)- The structure of the Sudakov exponent S{b,Q) is displayed in 
Eq. (IH). The spacelike distributions Vj ix,b) factorize at small b as 

"Pfi^^b) = [c;;f ®/J(x,6), (29) 

where C*"^^"^^ {x,b,as{bo/b)) are the spacelike Wilson coefficient functions for incoming par- 
tons. C^^ib) is a residual from the summation over the hadronic final states: 



CT\b) = 1 + ^^Cj. (-1 - yj + - ■ (30) 
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The right-hand side of Eq. (f27| depends on x via Pj {x,b). Nonperturbative contribu- 
tions to the z-Qow can be parametrized using the 6* prescription: 



W.{b, Q,x) = — Yl Cr*(6.) \CJJ, ® fa] {x, 6^)e-5p(^*,0)-5^p(^Q,-;^'.), (31) 

j=u,u,d,a,... 

where S^p{b,Q,x]b^) was found from an 0{as) fit [3| to the data 21, at 



8 X 10-5 < (x) < 0.11 and 10 < (Q^) < 2200 GeV^. The leading small-x behavior of 
Slip{b,Q,x;b^) is consistent with inverse power dependence: 

5^p(6,Q,a;;M ^ ^ + ... , (32) 



The best fit in Ref. [23| quotes Cq = 0.013 and p = 1. The x-dependence shown in Eq. (jT?^ 
describes, in particular, the data at g-r < 2 — 3 GeV, where the uncertainties due to the 
matching (switching) between the W and Y terms are small. By comparing two forms (f27|l 
and (|3T| for Wzip^Q^x), we find that the leading small-x term is a part of ix,b), in 
agreement with Eq. (fTTjl : 

Vf\x, b) ^ [c-f ) /.] (X, 6.) e-^-'^'(..M. (33) 

Theoretical uncertainties in the parametrization of V j (x, b) due to final-state fragmen- 
tation, matching, and higher-order corrections to the F-term are reduced by our choice 
of the experimental observable (2;-fiow) and kinematical region (g^ less than a few GeV). 
The exponential factor e-° °^^^^/^ parametrizes the full contribution of higher-order radia- 
five corrections of diverse wavelengths to V j (x, 6). It indicates that the SIDIS qx dis- 
tributions tend to become wider at smaller x, and it attributes the observed broaden- 
ing to scattering at b^ > 0.013/x. At x = 10~^ (10~^), this inequality translates into 
b^/b < 6o(a^/0.013)-^/^ ~ 1.3 (13) GeV, i.e., only scattering at low and intermediate energy 
scales is affected. Several intense mechanisms beyond 0{as) may turn on at such small x 
and 1/6 and cause the broadening. The rate may be increased at larger qx {1/b > 1 — 2 
GeV) by enhanced contributions from qq and qq hard scattering, as it happens in O^af) 
SIDIS cross sections at px = zqx ~ Q [23, 0, I^ll, Si ■ The rate may be reduced at smaller 
qx (l/&< 1 — 2 GeV) in the probed range of x because of the suppression of collinear split- 
tings of gluons by the logarithms ln(l/x) at higher orders of 2^. In addition to these 
leading-power (logarithmic in b) contributions, the x dependence may also arise from the 
unknown nonperturbative terms in Vj (x, b) proportional to the positive powers of b. In 
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the absence of a systematic computation of all such effects in the resummation formalism, 
the phenomenological broadening exponential serves as an approximate model of their total 
impact in the limited range 10"^ < a; < 10-2 covered by the HI data. The form ofV) '{x, b) 
is unknown at x outside of the HI data range, so that the parametrization (jHTni may not be 
extrapolated to x below 10"^. 



D. Relationship between resummation in DY and SIDIS 

The remaining step in the derivation of Eq. is the relationship (|THll between the 
timelike distributions Vj {x,b) and spacelike distributions Vj {x,b). This relationship 
follows from an elementary representation for spacelike (J = S) and timelike (/ = T) 



distributions 'P^'^^a;, 6) (42|: 



V';'{x,b)= nf\bo/b) U{bo/by^''Vj{x,b). (34) 

The function Hf\bo/b) = 1 + {as{bo/b)/7i)Hf'^^ + ... is composed of highly virtual cor- 
rections to the hard scattering vertex. The function ?7(6o/&)^^^ collects the soft subgraphs 
that are attached by gluon propagators to Hj^\bo/b). It is the same for / = S* and I = T. 
Vj{x,b) is related to the (n — 2)-dimensional Fourier-Bessel transform of the unintegrated 
(k-r-dependent) parton distribution function Pj(x,kr,C), 

Vj{x, b) = hm je-^'^^'^) j d"-2kre^^^ ''Pj"(x, kr, C)|. (35) 

Vj{x, kx, C) is given in terms of the quark field ipj by 

p,(x,kT,c) = EE / ^^7^^-'^'"''^'''^-'^m^ (36) 

spin color 

in a gauge n ■ A = with < 0. It depends on the gauge through the parameter 
( = (p ■ n)/\n'^\, where p'^ is the momentum of the parent hadron. 'Pj{x, b) is obtained from 
Vj{x, kx, C) by taking the limit C, oo. The partial Sudakov factor iS'(6, Q in Eq. (|35|l was 
derived in Ref. 0|. The function Vj{x, b) is the same in SIDIS and Drell-Yan processes Ql- 
Among the three terms on the right-hand side of Eq. (j34ll . only the hard vertex l-6p 
depends on the sign of the boson's virtuality. The functions and 1-6^ are given by the 
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same Feynman graphs in the different crossed channels, i.e., for bosons of timelike virtualities 
in the Drell-Yan case and spacelike virtualities in the SIDIS case. Therefore, 

nf\bo/b) 



Vf\x,b) 



nf\bo/b) 



vf\x,b)=r{b)VY\x,b) 



(37) 



where 



r{b) = 



nPibo/b) 



nf{bo/b) 



(38) 



A similar consideration for b = b* proves the relationship ([H^ between C"J^^\x,b^,) and 



^T/a. ^^^^ relationship (f2fl|) between c^'^\x,b;b^) and c*^'^^(x, 6; 6^,). 

According to Eq. (|3i|l . the dependence of the parton distributions Vj (x, b) and Vj {x, b) 
on X is determined by the same function Vj{x,b). More generally, the x dependence of the 
distributions V^^\x,b) {a = g,u,d, ...) is determined by the type of the parent hadron and 
parton flavor, but not by the type of the produced electroweak final state. For example, the 
X dependence of v'^P{x,b) is expected to be the same in the Drell-Yan process, 77 produc- 
tion via qq annihilation, and SIDIS; or in Higgs boson production and 77 production via gg 
fusion. This happens because all dependence of V^P [x, b) on the hard-scattering process is 
determined by the hard vertex function \'Ha\bo/b)\ (4^, which does not contain x. Va{x,b) 
are automatically process-independent in the representations for W{b,Q,XA,XB) that al- 



ways separate 



in 



from Va{x, b) (4 



. The best-fit phenomenological functions 



S^p'^^ {b, Q; b^) and S^p{b, Q; b^) were obtained in the CSS representation. For this reason, 
we will use this representation in the numerical analysis despite the additional complexities 
in the factorization relations. 



III. NUMERICAL RESULTS 
A. Overview 

In this section, we compare the resummed cross sections with the additional broad- 
ening term [p{x) 7^ 0] to the resummed cross sections without such a term {p{x) = 0]. We 
consider the decay modes eu, Z° — ee, and 77 and discuss the impact of the 

experimental acceptance cuts imposed on the decay particles. The resummation calculations 
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VS=14TeV VS=1.96TeV 




Figure 1: Parton momentum fractions xa and xb accessible in VF, Z, and Higgs boson production 
(for Mh = 120 GeV) in the Tevatron Run-2 {VS = 1.96 TeV) and at the LHC {VS = 14 TeV). The 
accessible ranges of xa and xb are shown by the solid lines. The contours of the constant rapidity 
y are shown by the inclined dotted lines. 

that include the decay of the vector and scalar bosons were described in Refs. and 0, 
respectively. The perturbative Sudakov factor was included up to 0{a1), and the functions 
pm{T) ^ 0{a) The numerical calculation was performed using the programs 

Legacy and ResBos P, 0|, and with the CTEQ6M parton distribution functions 3|- 

For the chosen parameters, the small-x broadening occurs when one or both longitudinal 
momentum fractions xa,b ~ Mye^'^ / \fS are of order, or less than, xq = 0.005. The accessi- 
ble ranges of xa, xb, and y at the Tevatron and LHC are displayed in Fig.[Tl Lower values of 
Xa can be reached at the price of pushing xb closer to unity, and vice versa. The number of 
events affected by the broadening can be determined from the rapidity distributions da/dy 
in Fig. [21'' In all scattering processes, most of the events occur at relatively small \y\. The 
rate in the forward regions is suppressed by the decreasing parton densities at a; — >^ 1. The 



The rapidity distributions and other observables inclusive in qr are insensitive to the qr broadening in 
our model, due to the Gaussian form of the broadening exponential assumed in Eq. (|2 
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Figure 2: Rapidity distributions da/dy for (a) boson production at the Tevatron and LHC; 
(b) boson production at the Tevatron and LHC; and (c) Higgs boson production at the LHC. 

broadening can be detected at all boson rapidities if > Mv/xq, or at forward rapidities 
\y\ > In Mv/{VSxo) , if the cross section is large enough to measure the Qt distribution 
at such \y\. The heavy bosons are identified in the experiment by observation of secondary 
particles from their decay. The magnitude of the broadening depends on the procedure 
applied to select the decay leptons or photons. The estimates for W bosons are modified 
by the integration over the unobserved rapidity of the neutrino, which mixes contributions 
from different ranges ofy. The impact of the boson's decay will be addressed quantitatively 
in the next subsections. 

In the Tevatron Run-2 (v^=1.96 TeV), the W and Z bosons are predominantly produced 
at xa ^ xb ~ My/VS > Xq. The broadening can only affect a relatively small fraction 
of the W and Z bosons with \y\ > In My /{^/Sxq) ~ 2, and it has negligible influence 
on most of the Tevatron observables. The feasibility of discovering qt broadening at the 
Tevatron depends on the available detector acceptance in the forward rapidity regions. The 
rapidity coverage of the D0 and CDF detectors is reviewed in the appendix. Our estimates 
suggest that verification of gr broadening at the Tevatron is viable in the near future. 

At the LHC (v^=14 TeV), at least one momentum fraction xa ov xb is of order, or less 
than, Xq in all W and Z boson events. Consequently the broadening is important at all 
rapidities. The magnitude of the broadening at the LHC generally differs between Z, W^, 
and W~ bosons, due to the differences in the bosons' masses and rapidity distributions. The 
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smaller mass of W bosons is conducive to the broadening. 70% of bosons are produced 
in ud annihilation, enhanced at \y\ > 1.5 by contributions from valence u quarks. Most of 
W~ and Z bosons are produced in scattering of valence d quarks and/or sea quarks, which 
tends to put more events in the central rapidity region. Consequently the fraction of forward 
rapidity events is larger in boson production than in W~ or Z boson production, and, 
of the three boson species, the bosons are the most sensitive to the broadening. The 
decay oiW^ bosons mixes contributions from different ranges of \y\, reducing the differences 
between qx broadening effects in the lepton-level observables. 

The Higgs bosons with the mass in the experimentally allowed range {Mh> 115 GeV) 
are produced via gluon fusion at somewhat larger momentum fractions, such as Mh/VS ~ 
0.0086 for Mh = 120 GeV. Their rapidity distribution in Fig.Efc) is narrower than in the W 
and Z boson cases due to the suppression of the forward rapidity regions by the decreasing 
gluon density g{x) at x 1. Consequently the broadening is reduced in the Higgs boson 
signal, but may be important in the background processes. 

Next, we present the quantitative results for the individual scattering processes. 

B. W and Z boson production in the Tevatron Run-2 

1. Z boson production 

The broadening may be most easily observed in the dilepton channel in Z boson pro- 
duction in the Tevatron Run-2. The strategy here is to exclude contributions from the 
central-rapidity Z bosons, which are almost insensitive to the broadening. If no distinction 
between the central and forward Z bosons is made (as, e.g., in the Run-1 analysis), the 
small-x broadening contributes at the level of the other uncertainties in the resummed form 
factor. Fig. infa) shows the Z boson distribution da/dqx, integrated over the Z boson ra- 
pidity y without selection cuts on the decay leptons. The cross section with the broadening 
term (dashed line) essentially coincides with the cross section without such a term (solid 
line). Both cross sections are dominated by contributions from x ~ Mz/VS ~ 0.046 ^ Xq, 
where the broadening function p(x) is negligible. 

In contrast, the small- a: broadening visibly modifies da/dqx at forward rapidities, where 
one of the initial-state partons carries a smaller momentum fraction than in the central 
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pp ^ Z^X ^ e^e X (Vs" = 1960 GeV) 
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Figure 3: Transverse momentum distributions of Z bosons in the Tevatron Run-2: (a) integrated 
over the full range of Z boson rapidity y; (b) for events with a large Z boson rapidity, | y | > 2; (c) for 
events with both decay electrons registered in the forward (ye+ > 2, yg- > 2, and Pre± > 20 GeV) 
or backward (ye+ < —2, y^- < —2, and PTe± > 20 GeV) detector regions. The solid curve is the 
standard CSS cross section, calculated using the BLNY parametrization of the nonperturbative 
Sudakov factor. The dashed curve includes additional terms responsible for the qx broadening in 
the small-x region, as in Eq. l(23|l . 
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Figure 4: Same as Fig. EIc), for the BLNY parametrization (solid), and various choices of the 
parameter xq in Eq. xq = 0.002 (dashed), 0.005 (dot-dashed), and 0.008 (dotted). 

region. Fig. IH^b) shows the cross section da/dqx for Z bosons satisfying \y\ > 2. The peak 
of the curve v^ith p{x) 7^ is lov^er and shifted toward higher qx- Fig.Efc) displays the cross 
sections with the acceptance cuts > 2, pTe± > 20 GeV or ye± < —2, pxe± > 20 GeV 
simultaneously imposed on both decay leptons. The cuts exclude contributions with \y\ < 2 
and retain a fairly large cross section 3.4 pb), most of which falls within the acceptance 
region (extending up to \ye\ ~ 3 in the D0 detector). The Run-2 can probably discriminate 
between the curves in Fig. [S^c) given the improved acceptance and increased luminosity of 
the upgraded Tevatron collider. 

To estimate sensitivity of the Tevatron Z boson production to uncertainties in the model 
for qt broadening, we vary the control parameter xq [cf. Eq. (f22|l ] within the limits compatible 
with the analyzed SIDIS z-Qow data. The parameter xq determines the upper boundary of 
the X range in which the broadening effects are important. As was observed in Ref. (2a|, 
the g-r-broadening was required in most of the rang e x G [10^^, 10~^] in SIDIS in order 
to obtain a good fit to d'Ez/d-QT from Refs. Given this observation, Xq = 0.005 

assumed in Fig. 01 is a conservative estimate of the upper boundary of the x region impacted 
by the broadening at HERA. Values of xq as high as 10^^ are compatible with the HERA 
data, while Xq ~ 10^'^ or smaller would lead to a poor fit to the HERA data in the range 

X e [xq, 10"^]. 
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If xo lies in the range 10~^ — 10~^, the broadening remains observable at the Tevatron. 
Fig. m shows the cross sections for xq = {0.002,0.005,0.008}, compared to the standard 
BLNY result. Although the broadening is reduced for xo = 0.002, it may still be discernible 
in a high-statistics data sample. On the other hand, the broadening can be also enhanced for 
Xo > 0.005, as demonstrated by the curve for xq = 0.008. The magnitude of the broadening 
for a given xq is affected by the parameter cq, which is anti-correlated with the power of x 
in the phenomenological factor e"'^"''^/^*' found in the SIDIS fit. Adjustments in cq can be 
compensated in some range by varying p and without deteriorating the overall quality of 
the SIDIS fit. For example, the central valu e cn = 0.013 was obtained for p = 1 in the best 



fit |2a| to both sets of the 2;-fiow data [23, |2^. The best fit to the earlier of the two HI 
data sets j^l was obtained in Ref. [3| for e~^NP = exp [—4.58 + 0.58/^x1, i-e., reduction 
in the negative power of x from 1 to 0.5 was partly compensated by the increase of Cq 
from 0.013 to 0.58. The magnitude of the broadening at the Tevatron for the superseded 
1/y/x parametrization is even larger than that shown in Fig. |4| Correlated uncertainties 
in the parameters of the SIDIS fit (and possibly a combined SIDIS/Drell-Yan fit) can be 
constrained in a future dedicated study, especially after improved SIDIS data from HERA-2 
become available. 



2. W boson production 

In the absence of leptonic cuts, W boson production is slightly more sensitive to the small- 
X broadening as compared to Z boson production because of the smaller mass and wider 
rapidity distribution of the W boson. The relevant qx distributions are displayed in Fig. [HI 
As in the case of the Z boson, the broadening is essentially negligible in the absence of selec- 
tion cuts on y or leptonic momenta [cf. Fig. Efa)]. In the forward region \y\ > 2 [Fig. E^b)], 
the broadening term shifts the peak of da/dgr significantly. The change 6qT ~ 750 MeV in 
the position of the peak is larger than the analogous change 6qT ~ 500 MeV in the Z boson 
cross section. 

Integration over the unobserved neutrino's rapidity in the eu decay channel mixes con- 
tributions from different rapidities of W bosons, thus enhancing the broadening at central 
electron rapidities and reducing the broadening at forward electron rapidities. The strongest 
effect occurs when only forward electrons are selected. Fig. ^c) shows the cross sections 
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Figure 5: Transverse momentum distributions for a combined sample of and W~ bosons in 
the Tevatron Run-2; (a) integrated over the full range of W boson rapidities; (b) integrated over 
the forward boson rapidities \ y\ > 2; (c) integrated with the selection cuts \ye\ > 2, pxe > 20 
GeV, > 20 GeV. The solid curve is the standard CSS cross section, calculated using the BLNY 



parametrization [2| of the nonperturbative Sudakov factor. The dashed curve includes additional 
terms responsible for the qx broadening in the small-x region, as in Eq. (j^Bl- 
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da/dqT with constraints imposed on the momentum of the electron (|t/e| > 2, pxe > 20 GeV) 
and missing transverse energy associated with the neutrino {]^t > 20 GeV). Although visi- 
ble, the difference between the solid and dashed curves in Fig. E^c) is less pronounced than 
in the W boson-level distribution in Fig. El^b) or in the 2^-boson lepton-level distribution in 
Fig. ini^c) . This is because some events in the sample with a forward electron rapidity and 
arbitrary neutrino rapidity are due to the decay of central-rapidity W bosons, which are 
not affected by the broadening. In contrast, events with two forward high-pT electrons in Z 
boson production are always due to the decay of forward-boosted Z bosons, and they are 
more sensitive to the broadening effects. 

Despite its modest magnitude, the forward qt broadening may affect the precision mea- 
surement of the W boson mass Mw in the Tevatron Run-2. The current goal of CDF and 
D0 experiments is to measure with the accuracy of 30-40 MeV per experiment. The 
W boson mass is commonly extracted from the distribution with respect to the leptonic 
transverse mass defined by 



in terms of the transverse momentum p^e of the electron and transverse momentum 
of the neutrino. In the experiment, the neutrino's transverse momentum is equated to 
the missing transverse energy, \ptu\ = $t- The determination of Mw from the shape of 
the kinematical (Jacobian) peak in da/dM^ is scarcely sensitive to the boson's transverse 
momentum g^-^ We have verified that the impact of the broadening on da /dM^ is negligible. 
However, reconstruction of I^t introduces a sizable uncertainty in the determination of M^- 
An alternative method is to extract My/ from the distribution da/dpxe- This method does 
not suffer from the complications associated with the reconstruction of J^t, but, as a trade- 
off, it is directly sensitive to the transverse momentum qt of the W bosons. In practice, the 
final uncertainties on Mw from the transverse mass and transverse momentum methods are 
comparable. As we will now show, the qt broadening can substantially affect Mw determined 
from the transverse momentum of electrons in the forward calorimeter. 

The distribution da/dpxe is shown in Fig. El Its Jacobian peak is located at 
Pxe ~ Mw /'^ ~ 40 GeV. To better visualize percent-level changes in da/dpTe associated with 

^ More precisely, a variation Sqr in the typical transverse momentum results in an error of order 
5M^^ /M^ ~ dq^/M"^ <C 1 in the typical transverse mass. 




(39) 
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Figure 6: Transverse momentum distribution da/dprp of the electrons from the decay of W bosons 
iji the Tevatron Run-2. CTEQ6M parton densities and BLNY nonperturbative Sudakov factor 
were used. 



in 



the broadening, we plot in Fig. [3 the fractional difference (da'^°'^ / dpTe) / [da'^^'^ / dpre) — 1 of 
the "modified" (mod) and "standard" (std) theory cross sections. The standard cross sections 
are obtained by taking p{x) = and some reference value of Mw, which we choose to be 
equal to 80.423 GeV. The modified cross sections are obtained by taking either p(x) 7^ 
(dotted curve) or a slightly varied W boson mass, Mw + 50 MeV (dashed curve) and Mw — 50 
MeV (dot-dashed curve) . The impact of qx broadening on the measurement of M]y can be 
judged by comparing the fractional differences due to the inclusion of the broadening term 
p{x) 7^ and explicit variation of Myy by ±50 MeV for p(x) = 0. 

A small increase of Mw results in a positive shift of the Jacobian peak, which reduces the 
cross section at pre < Mw/2 and increases the cross section at pTe > Mw A decrease in 
Mw shifts the Jacobian peak in the opposite direction. The broadening of da/dqx shifts the 
Jacobian peak in the positive direction. At |?/e| < 1 or in the rapidity-inclusive sample, the 
broadening can be mistaken for an increase of Mw by 10-20 MeV, as well as uncertainties in 
the parton distributions and/or nonperturbative Sudakov factor [cf. Fig. O^a)]. At \ye\ > 1, 
the small-x broadening exceeds the other theoretical uncertainties and is comparable with 
a variation of Mw by more than 50 MeV [cf. Fig. C^b)]. It remains observable even if the 
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Figure 7: The fractional difference in the distribution da/dpTe- (a) for the central-rapidity sample 
of electrons (|ye| < 1), and (b) for the forward-rapidity sample of electrons {\ye\ > 1). The shown 
curves are explained in the main text. 

electron's enerey is smeared due to the finite resolution of the detector. The experimental 
fractional resolution 6E/E near the Jacobian peak is typically less than 3% 0,0. We 
have estimated the energy smearing by convolving da/dpre with the Gaussian functions 
of the width corresponding to 6E/E between 1 and 5%. In all cases, the variation in the 
smeared da/dpxe due to the broadening has exceeded variations due to the shift of M^y by 
±50 MeV. 



C. Electroweak boson production at the Large Hadron Collider 

1. W and Z bosons 

At the LHC, the small- a; broadening may be observed in W and Z boson production 
at all rapidities. Very forward vector bosons probe x down to ~ 3 • 10~^, while the SIDIS 
data constrains the 6-dependent parton distributions at x above ~ 10~^. To stay in the 
kinematic region probed by the SIDIS data, we restrict our analysis to the events in the 
central-rapidity region, which we define by requiring the rapidities of the decay electrons to 
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Figure 8: Transverse momentum distributions of Z bosons at the Large Hadron Collider. The 
events are selected by requiring | ye I < 2.5 and "pTe > 25 GeV for both decay electrons. 

be less than 2.5. 

Fig. IHl demonstrates qt distributions of centrally produced Z bosons, selected by imposing 
the cuts luel < 2.5 andp^e > 25 GeV on both decay electrons. The distribution with p(x) 7^ 
is shifted toward higher qt- The broadening shift exceeds typical uncertainties in Snp- a 
similar shift would result from the increase of SNpip, Mz]K) from ^ 3.26^ in the BLNY 
parametrization to ~ 86^. The qx shift is even larger in production of W bosons [cf. Fig.lH], 
which we evaluate by requiring |t/e| < 2.5, pxe > 25 GeV, and > 25 GeV. The smaller 
mass Mw and less restrictive leptonic cuts result in smaller typical x probed by the W 
bosons, hence, in stronger broadening. The shift is slightly larger in boson production 
[cf. Fig. inia)] as compared in W~ boson production [cf. Fig. IHl^b)] because of the flatter 
y distribution for bosons [cf. Fig. Ela)]. The shown qt broadening propagates into 
the leptonic transverse mass and transverse momentum distributions. Both M^^ and pxi 
methods for the measurement of M]^ are affected in this case, in contrast to the Tevatron, 
where the M^'-' method is not susceptible to the broadening. 

2. Higgs boson 

The small- a; broadening is less spectacular, but visible, in the production of light Higgs 
bosons via the effective ggH vertex. Standard model Higgs bosons with mass below about 
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pp ^ W X ^ ev X (V^ = 14 TeV) 




2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 

[GeV] 

Figure 9: Transverse momentum distributions of (a) bosons and (b) W~ bosons at the 
Large Hadron Collider. The decay leptons are required to satisfy \ ye \ < 2.5, pTe > 25 GeV, 
> 25 GeV. 

140 GeV can be distinguished above the continuum background by their decay into a pair 
of highly energetic photons. The photon pairs produced from the Higgs boson decay have a 
larger typical qt than the photon pairs produced from the dominant background processes 
qq 77 and gg 77.^ For this reason, selection of photon pairs with g-r above 20-30 GeV 
can improve the statistical significance of the Higgs boson signal. 

The harder spectrum of the (^^f-dominated Higgs boson signal results from the larger 
leading-logarithm coefficient A^^'^ in the perturbative Sudakov factor (equal to Ca = 3 in 
gg channels and Cp = 4/3 = (4/9)Ca in qq channels). The Higgs boson distribution 

^ We refer to each hard subprocess according to its Born-level contribution, but the resummed cross sections 
include the relevant higher-order corrections to the Born-level diagram. For instance, the resummed cross 
section for qq 77 also includes the contributions from qq — > 771?, qg 779, etc. 
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Figure 10: The transverse momentum distribution da/dqT of on-shell Higgs bosons at the Large 
Hadron Colhder: (a) integrated over the full range of boson rapidities; (b) integrated over the 
forward region \y\ > 2. The meaning of the curves is explained in the text. 



peaks at = 10 — 20 GeV, as compared to ~ 5 — 10 GeV in gg-dominated Z boson 
distribution. A harder qx spectrum weakens the impact of the nonperturbative Sudakov 
function 0| and qt broadening. The broadening in the rapidity-inclusive Higgs boson cross 
section is further reduced because of the heavier Higgs boson mass and suppression of the 
forward- rapidity regions, as displayed in Fig. Efc). The nonperturbative and qr broadening 
terms in ^fgf-initiated processes are not constrained yet by the data. In the following, we 
evaluate the uncertainty in the Higgs boson g^ distribution by computing the resummed 
cross sections for several plausible choices of p{x) and SNp{b,Q;h^). 

Fig. ITOl a) shows g-r distributions for on-shell Higgs bosons with a mass Mh = 120 GeV 
and arbitrary rapidity. We first compare the cross sections for p(x) = and p(x) ^ 
(thick lines), with the same functions p(x) and SNp{b,Q; 6*) as in Z boson production. 
The difference between the two curves is minimal. The value of A^^^ refiects the dominant 
color fiow in the scattering. The bigger A^-^^ in Higgs boson production suggests that the Q- 
dependent part (and possibly other terms) in Snp niay be multiplied by a larger color factor 
than in the Drell-Yan process. To estimate this effect, we multiply Snp by Ca/Cf = 9/4. 
The resulting cross section (thin solid line) is close to the other cross sections. 
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The ln(l/x) terms in Higgs boson production may be enhanced as well, due to the direct 
coupling of the Higgs bosons to gluon ladders. However, the broadening would have to 
be quite large to affect qx of 20 GeV or more, i.e., near the location of the selection cuts 
on Qt of the photon pair. For example, increasing the function p{x) by a factor of two 
as compared to the Z boson case would lead to a distribution shown by the thin dashed 
line, which remains close to the other cross sections at qx above 20 GeV. The broadening is 
enhanced in the sample with \y\ > 2, which constitutes about 20% of the total cross section 
[Fig. fTOlf b)]. In this figure, the differences between the curves are more substantial at all qr- 
The theoretical uncertainties seen in Fig. irHTb) will have to be reduced to reliably predict 
Higgs boson production at forward rapidities. 

The small-x broadening will be also present in the main background channels, qq 
77 an d qq — » 77. Application of the CSS formalism in these channels was discussed in 



Refs. (S^, |22l, The dominant background process, qq 77, has the same structure of 
the resummed cross section as qq Z^. In analogy to Z boson production, the broadening 
effect in qq — > 77 is expected to be stronger than in gg or gg 77, so that the 

selection cut on qx of the 77 pair may have to be revised. Additional measurements will be 
needed to constrain the small-x behavior in the gg channel, e.g., by examining T production 
0, 3| or 77 production away from the Higgs signal region. As the mass of the photon pair 
decreases, the subprocess gg — > 77 becomes increasingly important, contributing about 40% 
of the total cross section at Q = 80 GeV at the LHC jzi|- By comparing qx distributions 
in pp — > 77 and pp ^ Z m the same region of Q, one may be able to separate the qq and 
gg components of the resummed cross section and learn about the x dependence in the gg 
channel. 



IV. CONCLUSION 



In this paper, we analyze consequences of possible non-trivial rapidity dependence in 
transverse momentum distributions of massive electroweak bosons at the Tevatron and LHC. 
A specific mechanism considered here assumes enhancement of radiative corrections in the 
impact-parameter-dependent parton distributions at x < 10^^. This enhancement can ex- 
plain the semi-inclusive DIS data from HERA, and it leads to wider transverse momentum 
distributions in the forward rapidity regions at hadron colliders. A measurement of the 
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rapidity dependence of qr distributions will test for the presence of such effects and verify 
universality of 6-dependent parton distributions, which was used in this paper to relate the 
SIDIS and Drell-Yan cross sections. 

Due to the two-scale nature of the transverse momentum resummation, large radiative 
corrections (such as non-resummed ln(l/a;) terms) may affect the Qt distribution even when 
they leave no discernible trace in the inclusive cross sections depending on one hard scale 
Q. The Collins-Soper-Sterman resummed form factor depends on the impact parameter b, 
which specifies the factorization scale ~ 1/b in the 6-dependent parton distributions. This 
scale changes from zero to infinity when b is integrated over the whole range in the Fourier- 
Bessel transform to obtain the resummed Qt distribution. When 1/b is much smaller than 
Q, the series ln"(l/x) in the CSS formula diverge at a larger value of x than the 

series a^{Q) ln"(l/x) in the inclusive cross sections. For this reason, the transition to the 
small-x dynamics may occur at larger x in qt distributions than in one-scale observables. 

For a realistic choice of parameters, the rapidity-dependent qr broadening may be discov- 
ered via the analysis of forward-produced Z bosons at the Tevatron Run-2. The broadening 
could be also searched for in the production of Drell-Yan pairs with a mass of a few GeV in 
fixed-target and collider experiments. Due to the different masses and shapes of the rapidity 
distributions for W^, and Higgs bosons, the magnitude of the broadening depends on 
the type and charge of the boson produced, and selection cuts imposed on the decay prod- 
ucts. In the Tevatron Run-2, the resulting change may shift the measured W boson mass 
in the pTe method by 10 — 20 MeV in the central region (|?/e| < 1) and more than 50 MeV 
in the forward region (|?/e| > !)• 

At the LHC, the most striking consequence is a harder qx distribution for W bosons 
displayed in Fig. IHl The predicted effect may easily exceed the other uncertainties in the 
resummed cross section, resulting in important implications for the measurement of the 
W boson mass from both transverse mass and transverse momentum distributions. The 
selection requirements imposed on the Higgs boson candidates in the 77 decay channel 
may have to be reconsidered to account for the non-uniform broadening in the signal and 
background processes. More generally, the broadening may bring about a revision of the 
earlier predictions that have neglected the correlation between the transverse momentum 
and rapidity in Drell-Yan-like processes. 
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APPENDIX: DETECTOR ACCEPTANCES AT THE TEVATRON AND LHC 

In this appendix, we review the angular coverage of the hadronic detectors and discuss 
feasibility of the observation of qx broadening at the Tevatron and LHC. 

In the Tevatron Run-2, adequate detector acceptance at large pseudorapidities 
Tj = — log(tan(^/2)) is crucial for the observation of qx broadening in Z boson produc- 
tion at y ~ |?7| > 2. Z bosons can be identified by registering two high-p^ electrons in the 
electromagnetic calorimeter, with simultaneous positive identification of at least one elec- 
tron by the tracking system. In the D0 detector, the silicon tracking system extends out 
to |?7| ~ 3 0,0,011 which provides an excellent opportunity for collecting a clean sample 
of forward-boosted Z bosons in events with two registered forward electrons. 

At CDF 0|, the unambiguous selection of the forward-boosted Z bosons is not possible 
due to the absence of tracking in the forward region. The CDF detector can efficiently 
track charged electrons in the main drift chamber at \ri\ < 1 and in the silicon system at 
|?7| < 2. Low efficiency tracking can be done up to |?7| ~ 2.6 by taking advantage of the 
longitudinally displaced collision vertex in a fraction of events. Energy deposits from the 
electrons can be registered in the plug calorimeter up to |?7| ~ 3.6. To increase sensitivity 
in the forward region, CDF could concentrate on events with one central (|?7| < 2) and 
one forward {\r]\ > 2) electron. Despite contamination by central Z bosons, this sample 
may reveal the small-x broadening once enough events are accumulated in the next years of 
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Run-2. 

The qt broadening in Z boson production can be also detected in muon decays. However, 
the coverage of the Tevatron muon systems (|?7| < 1.5 at CDF and \ri\ < 2 at D0 ) may be 
insufficient in the kinematic region affected by the broadening. 

Excellent opportunities to probe the x dependence of qr distributions will become avail- 
able at the LHC. Due to the larger center-of-mass energy, the qr broadening may affect W 
and Z boson production at all rapidities. In the ATLAS experiment, the tracking in the 
inner detector extends up to |?7| ~ 2.5. The central and forward electromagnetic calorime- 
ters cover \ri\ < 4.9, and the muon chambers cover |?7| < 2.7 In the CMS detector, the 
inner tracking system will register the high-py charged particles in the range |?7| < 2.5. The 
electromagnetic calorimeter extends up to |?7| = 3 Eventually, the TOTEM experiment 
will extend charged particle tracking and triggering capabilities of the CMS experiment into 
the range 3 < jr^l < 6.8 jisl. 



[1] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich, and G. Weiglein, Eur. Phys. J. C28, 133 
(2003). 

[2] R. K. Ellis and S. Veseh, Nucl. Phys. B511, 649 (1998). 

[3] A. Guffanti and G. E. Smye, JHEP 10, 025 (2000). 

[4] A. Kulesza and W. J. Stirling, Eur. Phys. J. C20, 349 (2001). 

[5] A. Kulesza, G. Sterman, and W. Vogelsang, Phys. Rev. D66, 014011 (2002). 

[6] X. Ji, J.-P. Ma, and F. Yuan, Phys. Rev. D71, 034005 (2005). 

[7] X. Ji, J.-P. Ma, and F. Yuan, Phys. Lett. B597, 299 (2004). 

[8] J. C. Collins, D. E. Soper, and G. Sterman, Nucl. Phys. B250, 199 (1985). 

[9] F. Landry, R. Brock, P. M. Nadolsky, and C.-P. Yuan, Phys. Rev. D67, 073016 (2003). 

[10] S. Berge, P. Nadolsky, F. Olness, and C.-P. Yuan (2004), hep-ph/0401128. 

[11] M. Dobbs et al. (2004), hep-ph/0403100. 

[12] Y. L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977). 

[13] V. N. Gribov and L. N. Lipatov, Yad. Fiz. 15, 781 (1972), [Sov. J. Nucl. Phys. 15 438 (1972)]. 

[14] V. N. Gribov and L. N. Lipatov, Yad. Fiz. 15, 1218 (1972), [Sov. J. Nucl. Phys. 15 675 (1972)]. 

[15] G. Altarelh and G. Parisi, Nucl. Phys. B126, 298 (1977). 



32 



[16] I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 

[17] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 44, 443 (1976). 

[18] E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977). 

[19] M. Ciafaloni, Nucl. Phys. B296, 49 (1988). 

[20] S. Catani, F. Fiorani, and G. Marchesini, Phys. Lett. B234, 339 (1990). 

[21] S. Catani, F. Fiorani, and G. Marchesini, Nucl. Phys. B336, 18 (1990). 

[22] G. Marchesini, Nucl. Phys. B445, 49 (1995). 

[23] J. C. Collins, Nucl. Phys. B396, 161 (1993). 

[24] R. Meng, F. L Olness, and D. E. Soper, Phys. Rev. D54, 1919 (1996). 

[25] P. Nadolsky, D. R. Stump, and C.-P. Yuan, Phys. Rev. D61, 014003 (2000), erratum: ibid., 
D64, 059903 (2001). 

[26] P. M. Nadolsky, D. R. Stump, and C.-P. Yuan, Phys. Rev. D64, 114011 (2001). 

[27] S. Aid et al. (HI Collaboration), Phys. Lett. B356, 118 (1995). 

[28] C. Adloffet al. (HI Collaboration), Eur. Phys. J. C12, 595 (2000). 

[29] B. A. Kniehl, G. Kramer, and M. Maniatis, Nucl. Phys. B711, 345 (2005). 

[30] P. Aurenche, R. Basu, M. Fontannaz, and R. M. Godbole, Eur. Phys. J. C34, 277 (2004). 

[31] M. Fontannaz, Eur. Phys. J. C38, 297 (2004). 

[32] A. Daleo, D. de Florian, and R. Sassot, Phys. Rev. D71, 034013 (2005). 

[33] M. Ciafaloni, D. Colferai, G. P. Salam, and A. M. Stasto, Phys. Lett. B587, 87 (2004). 

[34] J. C. Collins and A. Metz, Phys. Rev. Lett. 93, 252001 (2004). 

[35] J. R. Elhs, M. K. Gaillard, and D. V. Nanopoulos, Nucl. Phys. B106, 292 (1976). 

[36] M. A. Shifman, A. L Vainshtein, M. B. Voloshin, and V. L Zakharov, Yad. Fiz 30, 1368 (1979), 

[Sov. J. Nucl. Phys. 30 711 (1979)]. 

[37] S. Dawson, Nucl. Phys. B359, 283 (1991). 

[38] A. Djouadi, M. Spira, and P. M. Zerwas, Phys. Lett. B264, 440 (1991). 

[39] J. Kodaira and L. Trentadue, Phys. Lett. B112, 66 (1982). 

[40] C. T. H. Davies and W. J. Stirling, Nucl. Phys. B244, 337 (1984). 

[41] S. Catani, E. D'Emilio, and L. Trentadue, Phys. Lett. B211, 335 (1988). 

[42] R. P. Kauffman, Phys. Rev. D44, 1415 (1991). 

[43] C.-P. Yuan, Phys. Lett. B283, 395 (1992). 

[44] D. de Florian and M. Grazzini, Phys. Rev. Lett. 85, 4678 (2000). 

33 



[451 

I J 


S. 


Catani, D. de Florian, and M. Grazzini, Nucl. Phys. B596, 299 (2001). 


[461 

I J 


J. 


C. Collins and D. E. Soper, Nucl. Phys. B194, 445 (1982). 


[471 

I J 


J. 


C. Collins and D. E. Soper, Nucl. Phys. B197, 446 (1982). 


[481 

L J 


c. 


T. H. Davies, B. R. Webber, and W. J. Stirling, Nucl. Phys. B256, 413 (1985). 


[491 

I J 


G 


A. Ladinsky and C. P. Yuan, Phys. Rev. D50, 4239 (1994). 


[501 

I J 


R. 


K. Elhs, D. A. Ross, and S. Veseh, Nucl. Phys. B503, 309 (1997). 


[51] 


F. 


Landry, R. Brock, G. Ladinsky, and C.-P. Yuan, Phys. Rev. D63, 013004 (2001). 


[521 

I J 


F. 


J. Landry, Ph.D. thesis, Michigan State Universitv (2001), UMI-30-09135-mc (microfiche). 


[531 

I J 


J. 


Qiu and X. Zhang, Phys. Rev. D63, 114011 (2001). 


[541 

I J 


G 


P. Korchemsky and G. Sterman, Nucl. Phys. B555, 335 (1999). 


[55] 


S. 


Tafat, JEEP 05, 004 (2001). 


[56] 

I J 


R. D. Peccei and R. Ruckl, Phys. Lett. B84, 95 (1979). 


[57] 

L J 


R. 


D. Peccei and R. Ruckl, Nucl. Phys. B162, 125 (1980). 


[58] 

L J 


M 


. Dechantsreiter, F. Halzen, and D. M. Scott, Zeit. Phys. C8, 85 (1981). 


[59] 


J. C. Collins and D. E. Soper, Nucl. Phys. B193, 381 (1981), erratum: ibid., B213, 545 (1983). 


[60] 


C. 


Balazs and C.-P. Yuan, Phys. Rev. D56, 5558 (1997). 


[61] 

I J 


C. 


Balazs and C.-P. Yuan, Phys. Lett. B478, 192 (2000). 


[62] 


J. 


Pumplin et al., JEEP 07, 012 (2002). 


[63] 

I J 


J. 


Smith, W. L. van Neerven, and J. A. M. Vermaseren, Phys. Rev. Lett. 50, 1738 (1983). 


[64] 

I J 


T. 


Affolder et al. (CDF Collaboration), Phys. Rev. D64, 052001 (2001). 


[65] 

I J 


V. 


M. Abazov et al. (D0 Collaboration), Phys. Rev. D66, 012001 (2002). 


[66] 

I J 


C. 


Balazs, E. L. Berger, S. Mrenna, and C.-P. Yuan, Phys. Rev. D57, 6934 (1998). 


[67] 

I J 


c. 


Balazs, P. Nadolsky, C. Schmidt, and C.-P. Yuan, Phys. Lett. B489, 157 (2000). 


[68] 


p. 


M. Nadolsky and C. R. Schmidt, Phys. Lett. B558, 63 (2003). 


[69] 

I J 


A. 


Kulesza and W. J. Stirling, JEEP 12, 056 (2003). 


[70] 

I J 


E. L. Berger, J.-W. Qiu, and Y.-L. Wang, Phys. Rev. D71, 034007 (2005). 


[711 
['^1 


Z. 


Bern, L. J. Dixon, and C. Schmidt, Phys. Rev. D66, 074018 (2002). 


[72] 


S. Abadii et al. (D0 Collaboration), Nucl. lustrum. Meth. A338, 185 (1994). 


[73] 


S. 


Eagopian, in Proceedings of 7th International Conference on Advanced Technology and Par- 




tide Physics (Villa Olmo, Como, Italy, 15-19 October, 2001). 



[74] J. Elhson (D0 Collaboration) (2000), hep-ex/0101048. 



[75] F. Abe et al. (CDF Collaboration), Nucl. Instr. Meth. A271, 387 (1988). 

[76] ATLAS: Detector and physics performance technical design report, vol. 1 (CERN-LHCC-99- 

14) and vol. 2 (CERN-LHCC-99-15), 1999. 
[77] CMS electromagnetic calorimeter, Technical Design Report, CERN-LHCC-97-33 (1997). 
[78] V. Berardi et al., TOTEM: Technical design report, CERN-LHCC-2004-002 (2004). 



35 



